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Abstract. The influence of spacetime dynamics in hadronization via parton 

K^ \ coalescence at RHIC is investigated using covariant parton transport theory. 

/TT ' Key observables, the quark number scaling of elliptic flow and the enhancement 

^— V , of the p/tt ratio, show strong dynamical effects and differ from earlier results 

\^ • based on the simple coalescence formulas. 
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J . 1. Introduction 

^ , 

^ , Recent exciting puzzles in Au + Au reactions at y^J/vw = 130 and 200 GeV at 

RHIC are the quark number scaling of elliptic flow[n 121 HI and the weaker baryon 

suppression than that of mesons in the intermediate transverse momentum region 

^ i 2 < p_L < 5 GeV [H H H E). Parton coalescence [ El 0111111131111112113 is 

JH \ probably the most promising proposal at present to explain both phenomena. 

In the coalescence model, mesons form from a quark and antiquark, while 
baryons from three quarks or three antiquarks. The simplest version of the model is 
based on the "coalescence formula" (Eq. (1), or minor variation of it) that gives the 
hadron spectra in terms of the constituent phasespace distributions on a 3D space- 
time hypersurface. Remarkably, this simple approach can reproduce quite well the 
particle spectra at RHIC[|H1 EI and can explain the scaling of elliptic flow W2(pj.) 
with constituent number[ llU[ 111) . Furthermore, the parameterizations of the con- 
stituent distributions and the hypersurface are consistent with hadronization from 
a thermalized quark-antiquark plasma at RHIC. 

Nevertheless, these earlier studies left several important questions open. For 
example, it is known [ [TUl [TT] that the coalescence formula violates unitarity The 
yield in a given coalescence channel scales quadratically/cubically with constituent 
number, moreover, the same constituent contributes to several channels (including 
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fragmentation in certain schemes). In addition, it is not known whether the ex- 
tracted hadronization parameters (temperature, flow velocity, volume, hadroniza- 
tion time, etc.) are consistent with any dynamical scenario. In particular, the 
simple form of constituent phasespace distributions assumed ignores several kinds 
of phasespace correlations that would be present in a dynamical approach. 

The goal of this paper is to improve upon the above deficiencies and study 
how the dynamics of parton coalescence affects basic observables at RHIC, such 
as the p and tt suppression pattern {Raa) and elliptic flow V2{pi_) = {cos4>)pT- 
For simplicity, collisions at impact parameter 6 = 8 fm (« 30% centrality) are 
considered. Preliminary results were published in [I12|. 

2. Dynamical coalescence approach 

The parton coalescence formalism is largely based on studies of deuteron formation 
[ 1141 1151 1161 I17j . In nonrelativistic A^-body quantum mechanics, the number of 
deuterons (of given momentum) at time t is Ndit) — iVpairsTr[/5dy6(^^(i)], where 
Pd = \^d){^d\ is the deuteron density matrix, p'^'^\t) is the projection of the density 
matrix \^^^\t)){^^^\t)\ of the system onto the (two-particle) deuteron subspace, 
there are Npairs possible n—p pairs, while ^d and ^'-^'> (t) are the wave function of 
the deuteron and the system. The goal is a suitable approximation for the observed 
deuteron number Nd{t^ oo) because one cannot solve the full A^-body problem. 

One approach is to postulate that interactions cease suddenly at some time tf , 
and approximate p'-^'> in the Wigner representation as the product of classical phase 
space densities /„ x fp on the t = tf hypersurface. The approximation, at best, is 
valid for weak bound states and ignores genuine two-particle correlations. Applied 
to meson formation qq — > M, one obtains the simple coalescence formula 

^^^§^^9mJ n ('^'^'^'pO VFm(Ax, Ap) /,(pi, xi)/,(p2, X2) 5\p^p,-p2) (1) 

with Aa; = xi — X2 , Ap = pi ^ p2, and the meson Wigner function Wm [x, p) = 
J d'^b exp[—ibp]^*j^f{x~ b/2)^Mix + b/2). The degeneracy factor qm takes care of 
quantum numbers (flavor, spin, color). The formula for baryons involves a triple 
phasespace integral and the baryon Wigner function (which depends on two rela- 
tive coordinates and two relative momenta). The generalization to arbitrary 3D 
hadronization hypersurfaces is straightforward [ 1 1 51 117| . 

One difficulty with |^ is the proper choice of the hypersurface. In quantum 
mechanics the deuteron number is constant at any time after freezeout, however 
for free streaming ^ decreases [ITU [TC] with tf. Also, transport approaches (i.e., 
self-consistent freezeout) yield diffuse 4D freezeout distributions [ 1 1 81 I19L I2(J| . which 
cannot be well approximated with a hypersurface. These problems have been ad- 
dressed in [^J by Gyulassy, Frankel and Remler (GFR), where they derived a way 
to interface transport models and the coalescence formalism in the weak binding 
limit. 
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The GFR result is the same as ^, except that the weight Wm is evaluated using 
the freezeout coordinates (ti, xi), (t2, X2) of each constituent pair. When taking Aaf, 
the earlier particle needs to be propagated to the time of the later one, resulting 
in an extra term, e.g.. Ax = xi — X2 + (^2 — ti)vi if <i < ^2- The origin of this 
correction is that a weak bound state can only survive if none of its constituents 
have any further interactions. The generalization to baryons involves propagation 
to the latest of the three freezeout times. 

To investigate coalescence dynamics at RHIC, we implanted GFR into covari- 
ant parton transport theorv[ HTl 1221 1231 I19| . First the parton {g,u,d,s,u,d,s) 
evolution was computed until freezeout via the covariant Molnar's Parton Cascade 
(MPC) algorithm[ 23". For simplicity, only 2^2 processes were considered, with 
Debye-screened cross sections da/dt oc l/{t — fJ.%)^, fJ,D = 0.71 GeV. For the total 
gluon-gluon cross section two values were explored, agg — 3 nib (the typical pQCD 
estimate) and 10 mb. The quark-gluon and quark-quark cross sections were sup- 
pressed by the appropriate ratio of S'C/(3) Casimirs: Uqq — (4/9)crgg = (4/9)^cr£,g. 

At parton freezeout, the GFR formula was applied (in the two- or three-body 
center of mass frame) using, as common in transport approaches[ 1(3 , box Wigner 
functions W = H^j ^{^m - \xi - fj|)6(pm - \pi - p*,|), with Xm = 1 fm. (p„ is 
fixed by the normalization for W ^ This way (Q has a probabilistic interpretation: 
\iW = \ (and the quantum numbers match) the hadron is formed, otherwise it is 
not (W = 0). If several coalescence final states existed for a given constituent, one 
was chosen randomly with equal probability for all. Meson channels to tt, K , 77, ry', 
p, K* ^ w, $; and baryon channels to p, n, E, A, S, A, 51 were considered. Gluons 
were assumed to split to a q—q pair with asymmetric momentum fractions x ~ 
and 1, approximated as \g — + \q (or q). Easy color neutralization was also assumed. 
Partons that did not find a coalescence partner were fragmented independently via 
JETSET 7.4.10 [ 24,. JETSET was also used to decay unstable hadrons. 

Because Q requires a good sampling of full 6D phasespace, for the best statisti- 
cal accuracy, quarks and antiquarks were separately flavor-averaged ((u-|-(i-|-s)/3) at 
hadronization. This assumption of flavor equilibration, which is only approximate 
for strange quarks, will be lifted in the future (requires longer computation). 

The parton initial conditions for Au -(- Au at \/s = 200A GeV at RHIC with 
6 = 8 fm were taken from [ 22_. However, as illustrated for gluons in Fig. ^, for 
p± > 2 GeV LO pQCD minijet three-momentum distributions (solid line) were used 
(with a X-factor of 2, GRV98LO PDFs, and Q^=p|), which below pj_ < 2 GeV 
were smoothly extrapolated (circles) to yield a total parton dN{b~Q)/dy ~ 2000 
at midrapidity. This choice is motivated by the observed dNch/dy ^ 600 and the 
expectation that coalescence dominates the production. Perfect rj = y correlation 
was assumed. 

Figure QJi also shows that the initial condition reproduces fairly well the ob- 
served pion spectra in p + p at ^/s = 200 GeV at RHIC. Binary collision scaled 
ttq data from p + p [ I25j (crosses) are compared to the isospin-averaged pion spec- 
trum from the initial condition hadronized via either independent fragmentation in 
JETSET (dashed Une) or BKK95 fragmentation functions [EE] (dotted Hne). 
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a) b) 

Fig. 1. Results for Au + Au a.t yG= 200A GeV at RHIC with & == 8 fm. a) Initial 
gluon spectrum, and reproduction of pion spectra in p + p collisions at RHIC (see 
text for details); b) gluon quenching factor Raa as a function of pr for o'gg — 3 mb 
(open squares) and 10 mb (solid circles). 

3. Results on spectra and elliptic flow 

The final hadron momentum distributions are given by a convolution of three dy- 
namical effects, i) the evolution of parton phasespace distributions due to multiple 
scatterings, ii) the hadronization process, iii) resonance decays (hadronic final state 
interactions were ignored). 

The transport evolution quenches the parton spectra at high-p^ by about a 
factor of three {agg — 3 mb) to ten (10 mb), as shown in Fig. ^3- The ratio Raa 
of the parton spectrum at freezeout to the initial spectrum (which corresponds to 
no nuclear effects) is plotted. In the language of parton energy loss models [|22], 
the suppression in this study comes from incoherent elastic energy loss, while in the 
context of hydrodynamics, it reflects the cooling of the expanding system due to 
pdV work. These results are very similar to those in Ref. [|52] (see Fig. 9 therein), 
even though that study considered a pure gluon gas with thermal initial conditions. 

Figure El shows the influence of coalescence dynamics on the proton and pion 
nuclear suppression factor Raa, which is the ratio of the hadron spectra calculated 
from partons at freezeout to the spectra from the initial condition hadronized via 
independent fragmentation (same as binary-scaled p -\- p). If only fragmentation is 
considered, Raa for both species is about the same as that of partons (cf. Fig.^). 
On the other hand, parton coalescence enhances both pion and proton yields, and 
hence Raa, by as much as a factor of three in the "coalescence window" [E|EI| 1.5 < 
Pt < 4.5 GeV. The additional hadron yield comes dominantly from partons with 
0.5 < Pt < 2 GeV, as demonstrated in Fig. ^, where the fraction of partons that 
fragment independently (i.e., that find no coalescence partner) is plotted. About 
two-thirds of the partons, including essentially all partons above pt > 2.5 GeV, 
hadronize via fragmentation. 

Unfortunately, the inclusion of hadronization channels via coalescence does not 
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Fig. 2. Pion (left) and proton (right) nuclear suppression factor at RHIC in 
Au + Au at ^/s = 200A GeV, & = 8 fm, with parton coalescence (symbols) and 
without coalescence (curves with no symbols) , as a function of p±^ and parton cross 
section. Results for agg = 3 mb (dotted) and 10 mb (solid lines) are shown. 

solve the p/tt puzzle as evident from the striking similarity between the left and right 
panels in Fig. HI Though for mesons, the boundaries of the coalescence window agree 
well with estimates [ 8. 9: based on the simple coalescence formula JQ), for baryons 
the window is about the same as for mesons and does not extend to higher px- 
Furthermore, the enhancement is about the same for pions and protons (protons 
are only « 5% higher). Therefore the p/ii ratio at intermediate pT is the same as 
m. p + p, while the data shows an almost two- fold increase at this centrality[0]. In 
fact, the pion data[|2Hl favor (jgg k, 3 mb {R\j^ ~ 0.4), while the p/ttq systematics[ 
^ suggests R\a ~ 1 — 1-lj i-e-, Cgg « 10 mb. 

In Fig.|5Ji, the enhancement is characterized by Rcoah the ratio of the final spec- 
tra with hadronization via combined coalescence and fragmentation to that with 
hadronization via fragmentation only. Quite remarkably, Rcoai is almost indepen- 
dent of agg , despite the strong cross section dependence of the parton phasespace 
evolution shown by the quenching (cf. Fig. ^Jj) and elliptic flow[l22|. This very 
interesting aspect of parton freezeout would deserve a more detailed study. 

A qualitative explanation for the above features is that compared to 1^, in 
the weak-binding case assumed, baryon production is disfavored in a dynamical ap- 
proach. Baryons are more fragile than mesons because they have three constituents 
and therefore less chance to escape without further interactions. In other words, 
baryons are formed at later times on average, when the densities are smaller. This 
meson-baryon difference, which follows from the diffuse 4D nature of self-consistent 
decoupling[ ll8l I19L 1201 in spacetime, is absent if sudden freezeout on a 3D hyper- 
surface is postulated. To demonstrate this, we also plot in Figs. O^-b results for a 
(rather unrealistic) scenario with immediate freezeout on the formation r = O.lfm/c 
hypersurface, which does enhance the p/n ratio by a factor 1.5—1.7 and give a wider 
coalescence window for baryons. 

Figure 21 shows, for agg = 10 mb, the effect of coalescence dynamics on pion 
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Fig. 3. Results for Au + Au at ^s = 200 A GeV at RHIC with 6 = 8 fm for 
ffgg — 3 mb (dotted) and 10 mb (solid), or immediate freezeout at r = 0.1 fm/c 
(dashed-dotted line), a) Pion and proton enhancement from parton coalescence as 
a function of p± ; b) Fraction of partons that fragment independently as a function 
of PT- 



and proton elliptic flow, in particular on the scaling formula[|ni I10| 
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with constituent number n (n = 2 for mesons, 3 for baryons). The left panel shows 
that the V2 of direct pions and protons from dynamical coalescence is smaller than 
predictions based on ||2J), by 20 and 30%, respectively. The reduction is much larger 
than the few-percent corrections to ^ nonlinear in V2 or those due to higher-order 
flow anisotropies. It is a result of dynamical coordinate-momentum correlations 
that were ignored in earlier approaches that assumed x ~ p factorizable (or even 
spatially uniform) constituent phasespace distributions. 

The right panel in Fig. 01 shows that fragmentation contributions (and decays) 
further reduce the hadron V2. Unlike the amplification from coalescence, fragmenta- 
tion decreases (smears out) the anisotropy because hadrons from the parton shower 
are not exactly collinear with the originating parton (nonzero jet width, (1jj_1) > 0). 
The parton V2{p±) extracted from the hadron flows using (j^J would underpredict 
the real parton V2 by about a factor of two. Also, pions and protons differ on the 
scaling plot by at least 10%, maybe even as much as 30% (increased statistics to 
explore the pr/n > 1 — 1.5 GeV region is certainly desirable). A small 10 — 15% 
pion-proton flow scaling violation would be allowed by the published data [12 El- 
The high-statistics Run-4 data will hopefully provide much stronger constraints. 

Despite the flow amplification due to coalescence, the strong reduction of V2 
in the coalescence window, caused by the much smaller ~ 25 — 30% fragmentation 
yield, "threatens" to reopen the opacity puzzle at RHIC ['10'. The final magnitudes 
of proton and pion V2 are 30 — 50% below the data [^. Large parton cross sections 
(jgg w 20 — 30 mb, an order of magnitude above pQCD estimates, would be needed 
to generate a large enough anisotropy, at least in the approach presented here. 
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Fig. 4. Quark number scaled elliptic flow as a function of pT for pions (open 
squares) and protons (filled circles) in Au+Au at ^/s — 200 A GeV at RHIC with b = 
8 and agg = 10 mb, from hadronization via combined coalescence and fragmentation 
(right), and for primary hadrons (i.e., without decays) from coalescence (left). The 
constituent V2{pt) is also shown (solid lines). 

4. Conclusions 

The above findings demonstrate that coalescence is an important hadronization 
channel at RHIC. However, the results also show that dynamical effects on the 
baryon/meson ratios and elliptic flow scaling are potentially large. 

Further studies are needed to reveal what it takes to preserve the basic features 
of the simple coalescence formulas. It may be that the dynamics considered here was 
oversimplified, or the spacetime evolution in heavy-ion collisions is not understood 
well enough, or the QCD coalescence process cannot be approximated as a formation 
of weakly-bound states. 
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